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Staggered Chiral Perturbation Theory * 

C. Aubin a , C. Bernard a 

a Department of Physics, Washington University, St. Louis, MO 63130, USA 

We discuss how to formulate a staggered chiral perturbation theory (SXPT). This amounts to a generalization 
of the Lee-Sharpe Lagrangian to include more than one flavor (i.e., multiple staggered fields), which turns out 
to be nontrivial. One loop corrections to pion and kaon masses and decay constants are computed as examples 
in three cases: the quenched, partially quenched, and full (unquenched) case. The results for the one loop mass 
and decay constant corrections have already been presented in Ref. 



In order to reproduce the chiral behavior of 
staggered fermion simulation data as a 2 — » 0, 
one must account for the systematic effects aris- 
ing from 0(a 2 ) taste violations. (Taste refers to 
the staggered degrees of freedom resulting from 
doubling, while flavor refers to true quark fla- 
vor.) These taste violations are not negligible at 
current lattice spacings (a « 0.9 — 0.12 fm [5]). 
Lee and Sharpe have formulated [3] such a stag- 
gered chiral perturbation theory (SXPT) for one 
staggered field. Here we describe the generaliza- 
tion to the case of n flavors. 

We follow a three-step procedure: 

• First, we generalize the Lee-Sharpe Lagrangian 
to incorporate multiple flavors. This is the "4+ 
4+..." theory — n flavors with 4 tastes per flavor. 

• Next we calculate various meson properties; 
specifically we will calculate the Goldstone pion 
mass and decay constant for n = 3. 

• Finally, we adjust this result by hand to keep 
only one taste per flavor, thus accounting for 
taking the v / Det in simulations 

This can be done for any of the following cases: 
partially quenched (m va i cnco ^ 7n sca ), "full QCD" 
(Wvaienco = m sca ) and quenched (m valcncc ^ 
TO sea here we will show results for the 

partially quenched case. Complete results have 
been reported in Ref. 1 and preliminary fits to 
simulation data are shown in Ref. 0]. 

We begin here by formulating SXPT for n fla- 
vors, stating only the differences between the 
n = 1 H3 and n > 1 cases. We collect the Gold- 
stone bosons arising from the spontaneous break- 



down of SU(An) L x SU(An) R -> SU(An) vcc in the 
unitary field £ = exp(i$//), with <I> the hermi- 
tian An x An matrix (shown is the n — 3 case): 
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(1) 



where U = U a T a , K + = K+T a , etc. are each 4x4 
matrices, and the T a are the 16 taste matrices 
^5) £m> e ^ c - We include a An x An mass matrix 
Ad = diag(m u J, m^I, m s I, ...), with / the 4x4 
identity matrix. 

We have the Lagrangian 

C = f ^(d^d^) --nfTT(MH + MY^) 
o 4 



"•^singlet + CI V, 



(2) 



where the first two terms are the standard kinetic 
and mass terms from continuum XPT 5. £ s ingict, 
proportional to mj, arises from the anomaly, and 
gives the SU(An) singlets a large mass (we will 
take niQ — > oo in the end). The last term is the 
taste-symmetry breaking potential; for a single 
flavor, Lee and Sharpe found this to be the sum 
of six operators: 
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V = J2 C i° 

i=l 



tr(&E&E) + h.c] + ... . 



C7 2 l[tr(E 2 



(3) 
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The Oi are taste-breaking chiral operators that 
correspond to various four-quark operators. To 
get all terms in the potential in single-trace form, 
Lee and Sharpe performed a Fierz transformation 
on the original operators. 
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Because of this Fierz transformation, generaliz- 
ing to multiple flavors is tricky. Returning to the 
quark level, we recall that these operators come 
from four-quark operators with a net momentum 
change of 0(ir), which changes quark taste. This 
gluon exchange can also change color, but not 
flavor, so all four-quark operators must be of the 
flavor-unmixed form: q~i(j s ® £,t)qiq~j{js> <8> £t')?j> 
where (7., <g>£t) is the (spin<8>taste) matrix. In the 
naive theory, each bilinear is separately chirally 
invariant — these are the "odd" bilinears in the 
staggered theory. Thus, only the odd-odd four- 
quark operators are relevant here. 

Keeping all of the four-quark operators in the 
flavor-unmixed form, we then see that, using the 
standard spurion analysis, the taste matrices £t 
are singlets under the flavor SU (n) symmetry 
12] . This means we can make the replacement 

£t — •» (£t ) = £,t^ij, where i and j are flavor in- 

\ / ij 

dices. This must be done before the Fierz trans- 
formation performed by Lee and Sharpe to put 
the operators into single-trace form. Only then 
do the chiral operators follow from the flavor- 
unmixed four-quark operators. 

For the operators Oi, O3, O4 and Oq (which we 
combine into "W"), we can just replace £ t £t"^- 
Instead of the operators O2 and O5 we have 
four operators, O2V ,02A,Osv and O^a ( "W"), 
which are not in single-trace form. For exam- 
ple: 2 



±[Tr(^£)Tr(^ n) E) + /i.e.] and 



2V = J[ 

2A = ^(^^Tr^E) + h.c.]. The full po- 
tential is then V = U + W. 

One of the consequences of the two-trace form 
of the terms in hi ' is the appearance of quark-level 
hairpin terms similar to £ S mgict- These terms 
are of the form +-£{U t + D t + S t + ■ ■ -) 2 , where 
X' t depends on the taste channel we're discussing 
(t = V, A, or /): X'j = Ami/ 3 for the singlet 
case, and for the vector and axial-vector cases 
^v(A) = a "v{A)i w ith &'v(A) a nnear combina- 
tion of the coefficients in U 1 . 

These hairpins allow mixing among the flavor- 
neutral mesons (shown in Fig. [Q, which we can 
resum to give a non-diagonal propagator: 



G 



>MN 



X' 



AIN 



M 



(<1 2 



M 



m 2 N ) 



Figure 1. The U — U propagator at the chiral (a) and 
quark (b) level. Each x corresponds to an insertion 
of X' t , and the intermediate meson could be U, D or 
S (for n = 3). 
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(q 2 + m 2 K o)(q 2 + m 2 ) ■ 



(4) 



where the second term we denote by T>mn- 

It is interesting to note that the symmetries of 
the multiple-flavor theory are only slightly mod- 
ified from those of the single-flavor theory. For 
example: 

• The spontaneously broken symmetry (m = 
and a = 0): SU(4n) L xSU(4n) R -> SU(4n) vec , 
for any n > 1. 

• The residual chiral symmetry (m = and a 7^ 
0): U(n)i x U(n) r . We use £ and r (not L and 
R) to denote the left and right symmetries to 
remind us that they are mixtures of chiral spin 
and taste. 

• Fermion number [m 7^ (for n > 1, we as- 
sume all n masses are nondegenerate) and a 7^ 
0]: (7(1) vbc for one flavor, and this becomes 
(t/(l) V ec) n f° r more than one flavor. 

We show here the n = 3 partially quenched 
result for the pion mass and decay constant. 
Here a "pion" is any flavor-charged meson, and n 
refers to the number of sea quarks. We add two 
quenched valence quarks, x and y, and calculate 
the properties of the Goldstone "pion" P 5 + = xy. 
X and Y will refer to the flavor-neutral mesons 
composed of these quarks. 

For the pion mass, we calculate the pion self 



energy evaluated at p — 



All connected 



(at the quark level) terms cancel and we are left 
only with terms involving the disconnected terms. 
To one loop, we have: 



(m x + m y ) 



M<1 



1 



16tt 2 / 2 



d A q 



2V V XY 



2V 



XY 



-V 



XY 



a.t. 



(5) 



In this expression and below we leave off the ana- 
lytic terms (denoted by [a.t.]) for simplicity. Note 
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also that the pion mass vanishes in the chiral limit 
(m x , m y — y 0) as it must, since this is the true 
Goldstone boson. 

For the pion decay constant, we calculate the 
matrix element (0\j(i5\P$ (p)) = —ifp+P^, where 
jfjt, is the axial current for P 5 + . To one loop we 
find 



Sf P + 



n 1 

i 



(6) 



2V XY 



&>xx 



+8V% Y + <W\ Y ) + (V -> A) 



(7) 



The sum over Q is over mesons with one sea and 
one valence quark, and t runs over the 16 tastes. 

These expressions are quite general, and any 
relevant result can be found by taking limits be- 
fore performing the momentum integrals. For the 
"full QCD" case, set m x — m u and m y — rrid 
(m s ) for the true pion (kaon); for the quenched 
case, take m u ^ }S — » oo. 

Once one takes the desired limits, one must 
perform the integrals in Eqs. and Q. The 
integrands are ratios of products of terms of 



the form (g 



These can be expanded 



as sums of poles times their residues. Perform- 
ing the integrals, we keep only the chiral log- 
arithms (the analytic pieces are absorbed into 
"[a.t.]", and we leave off finite volume effects 
here). For a single pole, we get a term £(m 2 ) = 
m 2 ln(m 2 /A 2 ), while for a double pole, we get 
I(m 2 ) = - (ln(m 2 /A 2 ) + l) (with A the chiral 
scale). The residues multiplying these logarithms 
are complicated, and are given in full detail along 
with the analytic terms in Ref. pp. Also, see 
Ref. ^ for explicit expressions for m 2 K and f^. 

The last step before having the final result for 
the mass and decay constant is to adjust from 
four to one tastes per flavor. This is done using 
the "quark flow technique," in Ref. 1, where we 
determine where quark loops arise, and multiply 
each corresponding loop by 1/4. However, from 
that approach it is not clear that this works at 



all orders in perturbation theory. Further, even 
if it does work, the application of the quark flow 
technique could be quite complex. 

One can use the replica method to automate 
this: Take ni the quarks of each flavor i. Then 
calculate quantities of interest to a given number 
of loops as analytic functions of the rijj in the end, 
set each m = 1/4. This should take into account 
the transition from 4 — * 1 tastes per flavor at all 
orders automatically. 

Another interesting possibility that may arise 
in SXPT is that of an unusual phase. If (with 



a 2 A A = 



L s A 



'S : , 



$A < ^A. crit 



-4A, 



1 + a 2 A A /m| 5 

2 + 3a 2 A A /m 2 



(8) 



Sr, 



m 2 A could become negative for small, but non- 
zero m u = rrid- From fits, this does not appear 
likely for the physical case of QCD. A correspond- 
ing condition for m u — rrid = m s may be satisfied, 
although such a phase would disappear in the con- 
tinuum limit. The possibility of an unusual phase 
requires further study [Hj- 

Calculations for m 2 + and f p + are complete for 

the partially quenched, full and quenched cases. 
Preliminary fits are shown in Ref. [1], and ap- 
pear promising. The next goal is to include heavy 
quarks so as to calculate heavy-light decay con- 
stants, and an extension to baryons 
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